Novel Josephson Effect in Triplet Superconductor - Ferromagnet - Triplet 

Superconductor Junctions 
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We predict a novel type of Josephson effect to occur in triplet superconductor - ferromagnet - 
triplet superconductor Josephson junctions. We show that the Josephson current, /j, exhibits a rich 
dependence on the relative orientation between the ferromagnetic moment and the d- vectors of the 
superconductors. This dependence can be used to build several types of Josephson current switches. 
Moreover, we predict an unconventional temperature dependence of Ij in which Ij changes sign 
with increasing temperature. 

PACS numbers: 74.50.+r, 74.45.+C, 74.78.-w 



Josephson junctions made of unconventional supercon- 
ductors have attracted significant interest over the last 
few years due to their unconventional quantum trans- 
port properties 0, 0, 13, [3, HI • The latter are deter- 
mined by the formation of low-energy Andreev bound 
states I3, which, for example, lead to a low-temperature 
anomaly in the Josephson current 05 0- Josephson junc- 
tions consisting of triplet superconductors are of partic- 
ular interest since the superconductors' odd spatial sym- 
metry guarantees the formation of low-energy Andreev 
states for any orientation of the superconductors in the 
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FIG. 1: Schematic picture of a ID Triplet Superconductor 
- Ferromagnet - Triplet Superconductor (TSFTS) Josephson 
Junction. 

In this Letter, we predict a novel type of Josephson 
effect in one-dimensional (ID) triplet superconductor - 
ferromagnet - triplet superconductor (TSFTS) junctions 
(a schematic picture of such a junction is shown in Fig.^. 
In particular, we demonstrate a rich dependence of the 
Josephson current, /j, on the relative orientation be- 
tween the ferromagnetic moment, M, in the barrier and 
the dr^i^-vectors of the left and right triplet superconduc- 
tor ISfl, as described by the angles a and (see Fig. 
This dependence can be used to create Josephson cur- 
rent switches in which small changes of a or ^ can tune 
the junction between two "current states", in which Ij 



is either "on" (/j 7^ 0) or "off' (/j 0), or differs in its 
direction. These types of two-level systems are of great 
current interest in the field of quantum information tech- 
nology 10]. Moreover, we predict that a TSFTS junction 
leads to a qualitatively new temperature dependence of 
Ij such that in certain cases, Ij changes sign (i.e., its 
direction) with increasing temperature. Finally, we show 
that adiabatic changes in the orientation of d^^n yield 
a behavior of Ij which is significantly altered from its 
equilibrium form. 

We take the ID TSFTS junction to be aligned along 
the z-axis and be described by the Hamiltonian H = 
J dz dz' H{z^ z')^ where (Ti = 1) 
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and i^l{z) and '0cr(^) are the fermionic creation and an- 
nihilation operators for a particle with spin a at site z, 
respectively, a are the Pauli matrices, and A(z, z') = 
— A(z^ z) is the superconducting gap. The d-vector of 
the triplet superconductors on the left {j = z < 0) 
and right {j = z > 0) of the barrier is given by 
dj = (cos ^j, sin 6>j, 0). We take Ol = such that 
II X and Or = with d^^ lying in the spin xy- 
plane. The ferromagnetic junction, located at z = 0, 
possesses a moment Mq and represents a magnetic scat- 
tering potential, described by the last term in Eq.Q, 
with M{z,z') = Mo{cosa,sma,0)d{z)S{z'). 

We show below that two Andreev bound states with 
energies Ea,b are formed in the TSFTS junction. The 
Josephson current flows through these two states 11] and 
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In order to obtain the energies of the Andreev states, we 
start from Eq.Q and derive the Bogohubov-de Gennes 
(BdG) equations [HEi 

by introducing the unitary Bo- 

goUubov transformation 

'^t(^) = ^Un{z)an^Vn{z)bl^Wn{z)bn^Xn{z)al, 
n 

^i(^) = ^-Un{z)bn^Vn{z)al-Wn{z)an^Xn{z)bl, 
n 

where the sum runs over ah eigenstates of the junction. 
For the wave functions of the locahzed Andreev states on 
the left and right of the barrier, we make the ansatz 
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where Cj 



+1(— 1), j = L{R)^ and f<i~^ is the decay 
length of the Andreev state. Using for definiteness a 
p^-wave symmetry of the superconducting gap, A{kz) = 
A j sin kz (with lattice constant set to unity), the BdG 
equations take the form j{z) = j{z) where 
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In the following we take Al = Aq and Ar = Aqc'^^ for 
the superconducting gap on the left and right side of the 
junction. The above eigenvalue equation for E is subject 
to the boundary conditions ^l(O) = ^i?(0) and 



dz^R{0)-dz^L{0) = -2mMo 



P 
Pt 



^i?(0), (4) 



where P = as cos a — i&o sin a. The solution of the BdG 
equations yields two Andreev states with energies 



En. 



^ = kF^/DJA ^B-2C± 2y^{A-C){B-C), 



(5) 



Ao 
where 

A = cos2(0/2) [1 -Dsm'^{0/2)] , 

B = sin2((9/2) [l-L>cos2(0/2)] , 
C = [cos^((/)/2)-cos^((9/2)]cos^(a-6>/2), 

and D = [1 + ^^]"^ with g = mMo/kp. The two An- 
dreev states appear in the local density of states near the 




FIG. 2: (a) ±Ea,h and (b)/j as a function of 
^ = 0, D = 0.5 and several values of a. 



for T = 0, 



junction as two particle-like and two hole-like peaks at 
energies t£^o,6, respectively. 

We first consider di, || d^^, and present in Fig. [21 the 
energies of the Andreev states and the resulting Joseph- 
son current at T = as a function of (j) for several values 
of a. Note that for T = 0, only the negative energy 
branches of the bound states are populated, and thus 
contribute to Ij. The dependence of Ea^h on (j) is qual- 
itatively different for M || d^,/? and M ± ^l.r- For 
M ± dL,i? (o^ = 7r/2) the Andreev states are degener- 
ate with Ea^b/Ao = kp^D cos (j)/ 2^ and possess well de- 
fined spin quantum numbers ctq =t, (75 =|, since they are 
not coupled by the scattering at the ferromagnetic bar- 
rier. The zero energy level crossing at (j)Lc = (2n + l)7r 
(n=integer) occurs with dEa^h/d(j) 7^ 0, resulting in dis- 
continuous jumps of /j, as shown in Fig. [2Ib). It is in- 
teresting to note that Ea^h is identical to that of Andreev 
states near a potential scattering barrier In con- 
trast, for a 7^ 7r/2 the ferromagnetic barrier couples the 
Andreev states, yielding a splitting of their energies, as 
shown in Fig.[2Ia) for M || ^l^r = 0). This coupling 
yields dEa^b/dcj) = at (/)lc, such that Ij evolves con- 
tinuously with (j) [Fig. Eljb)]. While these results remain 
qualitatively unchanged with decreasing the depen- 
dence of Ij on D changes with the orientation of M and 
<^L,R' Specifically, for I) <C 1 one finds, to leading order 
in b, Ij - for a 7^ 0, but Ij - D^^^ for a = 0. 

Finally, note that for = and any a, the induced An- 
dreev states with Ea^b/ An = kp^D are identical to the 
impurity (Shiba) states 14] induced by a single (static) 
magnetic impurity in a ID triplet superconductor. 

The Josephson current in a TSFTS junction exhibits 
an unconventional temperature dependence in that Ij 
can change sign, and thus its direction, with increasing 
temperature, as shown in Fig.jSfa) for ^ = and (/) = 7r/2 
(we assumed a BCS temperature dependence of the su- 
perconducting gap). In order to understand this sign 
change, we consider the (/)- dependence of Ea^b shown in 
Fig.[2fb) and note that at T = 0, only the branches indi- 
cated by 1 and 2, belonging to Andreev states a and 6, re- 
spectively, are occupied. Since the derivatives dEa^b/d(j) 
possess opposite signs for the two Andreev states, the cor- 
responding currents through them, /j < and /j > 0, 
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FIG. 3: (a) Ij as a function of T/Tc for 6* = 0, = O.Stt and 
D — 0.2. (b) d:Ea,b as a function of (j) for a = 0. The dotted 
line corresponds to the parameters in (a). 



flow in opposite directions with Ij > With in- 

creasing temperature, the occupation of branches 2 and 3 
changes more rapidly than those of branches 1 and 4. As 
a result, the magnitude of Ij decreases more quickly than 
that of /j, and the total current, /j, eventually changes 
sign. Moreover, the qualitative nature of the temperature 
dependence can be altered via a rotation of M. Specif- 
ically, as a increases, the T = value of Ij decreases 
while Ij remains practically unchanged. This leads to 
a qualitative change in the temperature dependence of 
Ij such that for a > 0.2l7r, Ij does not undergo a sign 
change with increasing temperature [Fig. Olja)] . 

The unique dependence of Ij on (j) and on the orien- 
tation of M can be used to build a Josephson current 
switch in which Ij is turned "on" or "off' via the rota- 
tion of M. This is shown in Fig.^a) where we present Ij 
as a function of a for di, || di^ at T = 0. Note that the 
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FIG. 4: (a) Ij and (b) ±Ea^b, as a function of a for D 
at T = 0, ^ = and several values of 6. 
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Qf-dependence of Ij changes signiflcantly with increas- 
ing 0. In particular, Ij becomes sharply peaked around 
a = (2n + l)7r/2 with integer n (i.e., for M ± ^l,r) 
as (j) approaches tt. As a result, small variations in a 
lead to large changes in the magnitude of the Joseph- 
son current, and can thus tune the junction between an 
"on" -state (/j 7^ 0) and an "off' -state (/j ^ 0). This 
behavior of Ij is directly reflected in the a-dependence 
of the bound state energies, as shown in Fig.^b). While 
for = the Andreev states are a-independent (see 
above), Ea^b oscillate sinusoidally with a for (j) = 7r/4, 
as shown in Fig. ^b). As approaches tt (see, e.g., 
^ = O.Qtt), the energies of the Andreev states alternate 



in being close to zero and almost a-independent for half 
of the period and sinusoidal for the other half. Finally, 
for ^ = TT (not shown) the bound state energies are given 
by Ea/ Ao = 2kF\/D(l- D) cos a and E^ = 0. Note 
that for (/) < TT, the Andreev states exhibit no zero-energy 
crossings when M is rotated . 

Another type of Josephson switch can be created by 
rotating dn for flxed M and di,, as shown in Fig. [3 for 
T = 0. For a = 0, ^ = 7r/2 and D = 0.7, Ij exhibits 
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FIG. 5: Ij as a function of 6* for T = 0, D = 0.7, several 
values of 0, and (a) a = 0, (b) a = 7r/4, and (c) a = 7r/2. (d) 
Ij for D 0.2 and a = 0. 

an almost perfect square wave form, is symmetric around 
= niT^ and is nearly ^-independent between the discon- 
tinuous jumps, as shown in Fig. (SJa). For cj) 7^ 7r/2, this 
symmetry is broken, and the ranges of over which Ij is 
positive or negative are unequal, but the discontinuous 
jumps in Ij persist. For a = 7r/4, the Josephson current 
is skewed, and in some regions varies nearly linearly with 
[Fig. ISIb)]. For a = 7r/2 and (j) = 7r/2, Ij is again 
symmetric around = utt [Fig. El^c)]. With decreas- 
ing the Josephson current becomes 6>-dependent be- 
tween discontinuous jumps, even for a = and = 7r/2 
[Fig-Efd)]. Thus small changes in the relative alignment 
of the d- vectors can switch the junction between two cur- 
rent states with opposite direction of Ij. 

A third type of Josephson switch can be formed via a 
correlated rotation of M and d^^. Speciflcally, we note 
that Ea^b^ Eq.(|S|), depends on a only via cos^(a — 0/2). 
We therefore consider a simultaneous rotation of d^^ and 
M such that a = ^/2, and present the resulting Ij as a 
function of in Fig. Efa). Ij then exhibits square wave 
oscillations with discontinuous jumps from /j to a 
negative (positive) value for < ^ < tt (tt < < 27r). 
To understand these sharp transitions, we plot the (j)- 
dependence of —Ea,b in Fig. Efb) (only these branches 
contribute to /j at T = 0). For = tt the Andreev states 
are degenerate, with dEa^b/dcj) < 0(> 0) for < (/) < tt 
{tt < (j) < 27r). In contrast, for = 7r/2 the degeneracy 
of the Andreev states is lifted with dEa/dcj) ^ —dEb/d(j) 
in some range of (j) (shown as dashed lines in Fig. Efb)) 
such that /j ~ 0, but /j 7^ in other regions of (j) (shown 
as solid lines). Note that as (j) increases from zero, the 
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FIG. 6: Correlated rotation of M and with a = 0/2 and 
D = 0.01 at T = 0. (a) Ij as a function of 0, and (b) -Ea,b 
as a function of 0. 

range in in which Ij increases, while the magnitude 
of I J decreases. This correlated rotation thus represents 
another type of Josephson current switch in which Ij can 
be turned "on" or "off' . 

Above, we computed Ij assuming that the phase differ- 
ence between the superconductors, ^, or the relative ori- 
entations between M and dL,R (i.e., a and 0) are main- 
tained for a sufficiently long time such that the system is 
in thermodynamic equilibrium. The qualitative form of 
I J changes, however, if M or d-L^n are rotated adiabati- 
cally on a time scale which is shorter than the relaxation 
time required for the occupation of the Andreev states 
to achieve their equilibrium values. In this case, the oc- 
cupation of the Andreev states remains unchanged from 
an initial value during the adiabatic rotation 15]. The 
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FIG. 7: (a) /j at T = for an adiabatic rotation of M and 
in equilibrium (see also Fig. |5Ja)) as a function of ^. a = 0, 
= 7r/2 and D = 0.7 . (b) ±Ea,b as a function of 0. The 
black solid lines represent the occupied states if is changed 
adiabatically. 

Josephson current at T = for an adiabatic rotation of 
di^ and in the equilibrium case are shown in Fig.[71^a). To 
understand the qualitative form of /j, we note that in the 
equilibrium case, only the negative energy branches of the 
Andreev states contribute to Ij. In contrast, if is adi- 
abatically changed (from the state with ^ = that is in 
thermodynamic equilibrium), then the energy branches 
of the Andreev states represented by the black solid lines 
in Fig. [7fb) remain occupied, while those shown as red 
dashed lines remain unoccupied. As a result, Ij in the 
adiabatic case changes continuously with and does not 
exhibit the discontinuous jumps shown by the equilib- 
rium Ij. Note that Ij is 27r-periodic in both cases. A 



special case arises for a = 7r/2, when the periodicity of 
I J is changed to An (not shown). An adiabatic change of 
(j) also leads to a 47r-periodicity of /j, in analogy to the 
fractional ac- Josephson effect discussed in Ref. p . 

Finally, scattering off the ferromagnetic barrier in gen- 
eral leads to a suppression of the superconducting order 
parameter near the barrier, which was not accounted for 
in the above approach. However, in nodeless supercon- 
ductors, such as the one discussed above, the order pa- 
rameter recovers its bulk value on a length scale set by 
which is assumed to be much shorter than the 
decay length of the Andreev bound state 0, • As a 
result, the spatial variation of the order parameter leads 
to only weak quantitative and no qualitative changes in 
the induced fermionic bound states 16]. We thus expect 
that the results presented above are unaffected by the 
inclusion of the order parameter suppression. 

In summary, we predict a new type of Josephson effect 
in TSFTS junctions, in which Ij exhibits a rich depen- 
dence on the relative orientation of M and d^.^^^. This 
dependence can be used to build several types of Joseph- 
son current switches in which small changes of a or ^ can 
tune the junction between two current states. We predict 
an unconventional temperature dependence of Ij such 
that for certain orientations of M and d^,/?, Ij changes 
sign with increasing temperature. 
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